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In fractional quantum Hall physics, the Hilbert space is projected to a single Landau level and the 
entire Hamiltonian consists of just the projected inter-electron interaction. Haldane's pseudopoten- 
tial formalism has been an extremely useful tool both for understanding these interactions, and for 
understanding the quantum Ifall states that result. In the current paper we consider the analogue 
of this pseudopotential construction that results from general Af-body interactions rather than the 
usual (Coulomb) two-body interaction. 



I. INTRODUCTION 

In high magnetic fields, the electronic states of two di- 
mensional systems break up into degenerate bands known 
as Landau levels^. In high enough fields, i.e., in the frac- 
tional quantum Hall regime, one assumes that the spac- 
ing between these levels is large enough such that all of 
the degrees of freedom are restricted to a single Lan- 
dau level. In this case, the entire Hamiltonian of the 
system is simply the projected inter-particle interaction. 
In an extremely insightful early paper by Haldane^i^, it 
was pointed out that any translationally and rotation- 
ally invariant two-body interaction, projected to a single 
Landau level, could be completely described by a set of 
"pseudopotential coefRcients" (Vl with L > 0) describ- 
ing the energy cost to have a pair of particles in a state 
of relative angular momentum L (for electrons L must be 
odd whereas if we were considering quantum Hall effect 
of bosons^ L would have to be even). This formalism 
turned out to be useful not just for describing the details 
of the interaction, but also for describing the resulting 
fractional quantum Hall states. For example, the Laugh- 
lin V = 1/3 state is precisely described as the unique 
highest density zero energy state of a potential which has 
Vi positive and Vl = for all L > 1. Furthermore, using 
the pseudopotential formalism, an interacting system in 
a partially filled higher Landau level can be mapped to 
an equivalent problem in the lowest Landau level with a 
modified interaction^. 

In recent years, the study of the fractional quantum 
Hall effect has begun to consider interactions beyond 
just Coulomb two-body interaction. Indeed, a particu- 
larly interesting class of quantum Hall states, including 
the Moore-Read state- and the Read-Rezayi^ states are 
exact highest density zero energy states of M-body inter- 
actions with M > 2. More recent work^i^has considered 
even more complicated Af-body interactions. These more 
complicated many particle interactions may not just be 
theoretical tools. For example, it has been shown^ how 
multi-particle interactions may be engineered in trapped 



atom systems. More generally, multi-particle interac- 
tions could easily result from integrating out virtual ex- 
citations of any sort. Furthermore, while the above 
mentioned quantum Hall states are exact ground states 
for very specific many-particle interactions, they may 
also be extremely good trial states for more realistic 
interactions^^. More importantly, in many cases, these 
exact ground states provide a particularly simple (and 
easily studied) representative of an entire phase of mat- 
ter. 

Despite the increasing interest in these multi-particle 
interactions, generalizations of Haldane's formalism have 
not, to our knowledge, been systematically constructed. 
(Note, however, that some progress in this direction has 
been given in Ref . [ill in a somewhat different language) . 
Performing this construction explicitly is the main objec- 
tive of the current paper. 

The outline of this paper is as follows. In section |TT] 
we start with a brief review of Haldane's pseudopoten- 
tial construction for simple two-body interactions. In 
Section IIIII we consider how this construction should be 
generalized to Af-body interactions. We find that the 
construction is more complicated than the two-body case 
because the simple pseudopotential coefficients need to 
be generalized to pseudopotential matrices. Section HVl is 
devoted to calculating the dimension of the resulting ma- 
trices, and specifying a convenient basis in which to ex- 
press these interaction matrices. In section IVl we briefly 
discuss how these results might be used and why they 
are interesting. For clarity, some simple examples of us- 
ing these results for three-body interactions are given in 
detail in the final appendix. 



II. REVIEW OF HALDANE'S CONSTRUCTION 

Before launching into the general construction for Af- 
body interactions, it is useful to review how Haldane's 
formalism works for two-body interactions^*^ (albeit in a 
slightly different language than that of Haldane's original 
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construction). Projected to a single Landau level, the 
Hamiltonian consists of an interaction term only 

H = J2^{r,-r,) (1) 

It will be crucial that the interaction term is both rota- 
tionally and translationally invariant. 

Since the Hamiltonian operates on two electrons at a 
time, we will consider two electron wavefunctions of the 
form \E'(ri,r2). This is sufficient to analyze the effect 
of the Hamiltonian in general since any many electron 
wavefunction can be decomposed into a sum of a two 
particle wavefunction times a wavefunction of everything 
else. For example, to find the effect of the term F(ri — 
on a general multiparticlc wavefunction we can always 
decompose 

^'(ri,...,rw) -Xl*''(''i'''2)^a(r3,...,r7v) (2) 

a 

and clearly the V{ri — term of the Hamiltonian only 
operates on the term and not the 4'a term. Thus, we 
need only focus on two electron wavefunctions ^(r^, rj). 

It is then useful to further decompose the two electron 
wavefunction into center of mass and relative degrees of 
freedom 

*(r„r,) = ^A.vpfM f'-l±Il\ (r, - r,) (3) 

b,c ^ ^ 

where vj/p^f and 4"™' form a basis for the center of mass 
and relative wavefunctions respectively. We note that the 
Hamiltonian, which is translationally invariant, acts only 
on the relative coordinate. 

In the lowest Landau level, using the usual analytic 
notation for wavefunctionsi where single particle orbitals 
are given by z™ exp(— |zp/4), and the magnetic length is 
taken to be unity, we can thus rewrite any product of the 
wavefunctions 5*^*^ and ^'JJ^' as 

where we have separated out the Gaussian factors. Both 
ip'-'-'^ and ijj^'^'' must be analytic polynomials, and some- 
times we will notate ip^^'' as simply ip. Throughout this 
paper we will be discussing the planar (disk) geometry. 
In Appendix |A] other geometries (sphere and torus) are 
briefly discussed. 

A convenient complete basis to describe the possible 
relative wavefunctions is now given by 

|L;i,j) = Cl(z,-z,)^ (5) 

where L is the relative angular momentum between par- 
ticles i and j and Cl is an appropriate normalization 
constant. Note that if the particles are fermions then L 



must be odd, and if they are bosons, then L must be even. 
We can now insert this complete set into the Hamiltonian 

^ = EEI^'*'J') {L;hj\Vir,-Vj)\L';t,j) 

i<j L.L' 

(6) 

Note that since the interaction is translationally invari- 
ant, we need only insert a complete set for the relative 
degrees of freedom (I.e., the interaction is diagonal and 
trivial in the center of mass degree of freedom) . 

The rotational invariance of the interaction now makes 
the matrix element diagonal in L so we obtain the Hal- 
dane Hamiltonian^i^ 

where 

P,'^^=\L:t,j){L;t,j\ (8) 

is a projection operator that projects particles i and j 
to a state of relative angular momentum L (within the 
lowest Landau level) . In Eq. [7l 

Vl,2 = {L;i,j\V{r,~r,)\L-i,j) (9) 
- {L\V\L) (10) 

is known as the pseudopotential coefficient, and we have 
added the subscript 2 here to indicate that this is a two- 
body interaction. In Eq. [TO] we have written this matrix 
element in a convenient shorthand, since by symmetry 
between particles, Vl,2 is independent of which i and j 
is chosen. 

One of the great advantages of the pseudopoten- 
tial formalism is how easily it generalizes to higher 
Landau levels^'"'^?' as well as to even more complicated 
situationsi^. The simple generalization stems from the 
fact that in any Landau level, there is a one-to-one map- 
ping to another system completely in the lowest Landau 
level with a modified inter-electron interactio n^'^^i^^ . To 
see how this mapping works, we consider a system in the 
first excited Landau level (ILL). Here, one can write any 
single electron state as a Landau level raising operator 
applied to a lowest Landau level state. In this way, 
an orthogonal set of relative angular momentum states 
analogous to Eq. [5] may be written in the ILL 

\L;t,j-lLL)^ala]\L;t,j) (11) 

wherein 

al = V2(-a.. -I- z*/A) (12) 

where the derivative acts on the Gaussian factors as well 
as on the polynomial part of the wavefunction (and * 
means complex conjugation). Using this orthogonal set 
for states in the ILL, and following the same argument, 
we obtain exactly the same Haldane Hamiltonian (Eq. [T]) 
but with modified pseudopotential coefficients 

Vl,2 = {L;i,j\aiaj V{r, - rj) a|a]|L;i,j). (13) 
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Further discussion of Landau level raising (for the case 
of three-body interactions) is given in example 4 of ap- 
pendix |e1 

III. M-BODY INTERACTIONS 

We now generalize Haldane's argument to M-body in- 
teractions. We thus consider a general M-body interac- 
tion Hamiltonian 

H= V{v^,...,vm) (14) 

il<i2<---<iM 

We note that, although not a requirement, we will typi- 
cally want to think about Af-body interactions that are 
not reducible to sums of interactions of a smaller num- 
ber of particles. This issue is discussed in more depth in 
appendix |B] 

Again crucial to our construction is the assumption 
that V is both translationally and rotationally invari- 
ant. We then consider a general M-body wavefunction 
^'(ri, . . . , tm), which we rewrite in terms of relative and 
center of mass coordinates 

«'(ri,...,rM) = ^Ab,c (15) 

-I^r( ^^+^^^; + ^" )^r(rr,...,rM) 

The statement that 5''"'^' is a wavefunction for relative 
motion means it is translationally invariant, or 

vI/'-^'(ri,r2,...,rM)= (16) 
^''''='(ri+x,r2+x...,rM+x) 

for arbitrary x. 

Once again, we can assume we are in the lowest Landau 
level so we can work with analytic forms of the wavefunc- 
tion, so we write the product \I/'^*^v[/'''^' as 

V'^"' (ir E -li ^) ^(-1, . . • , -M)e-i i-i^ (17) 

where is then the relative wavefunction, which must 
be a translationally invariant polynomial (and must be 
antisymmetric for fermions and symmetric for bosons). 
As in the two-body case, we can categorize these rel- 
ative wavefunctions in terms of an angular momentum 
quantum number. We define the "relative angular mo- 
mentum" operator on the relative wavefunction ip to be 
the total degree of the complex polynomial V'- If the 
polynomial ?/; is homogeneous of degree L, it is a rel- 
ative angular momentum eigenstate with eigenvalue L. 
Once again we note that our construction here pertains 
to the planar (disk) geometry. The sphere and torus are 
discussed briefly in appendix \K\ 

As in the two-body case, we are now seeking a complete 
basis in which to write our relative wavefunctions. We 



can certainly write any relative wavefunction ip as, a, sum 
of components with all possible values of L > 0. Le., 

■4){zi, zm) = '^tphizi, ■ ■ .,zm) (18) 

L 

where i/'l has relative angular momentum L, i.e, is ho- 
mogeneous of degree L. Unfortunately, as compared to 
the two-body case, it is no longer true that for a given L 
there is only a single possible wavefunction. Indeed, as 
we will see in great detail in section ITVl below, there may 
be many M-body wavefunctions with the same relative 
angular momenta. Thus, within the space of states with 
relative angular momentum L, we arbitrarily define an 
orthonormal basis which we write as \L^q). The dimen- 
sion d of this space (and hence the number of different 
values of q) depends on L, M and whether we are consid- 
ering bosons or fermions. In section lTVl we will determine 
these dependencies. 

As in the two-body case, since the interaction is inde- 
pendent of the center of mass of the M particles, we can 
then insert this complete set of relative wavefunctions 
into the Hamiltonian to obtain 

^ |L,g;ii, . . . ,Uf)V^;'^(L,9';ii, . . . 

il<...<iM L.q.q' 

(19) 

which is a sum over all possible M-clusters. Here the 
matrix element is a Hermitian matrix 

^L,M = {L,q]h, ■ ■ ■ ,iM\V{vi, . . . ,YM)\L,q']ii, . . . ,iM) 
= {L,q\V\L,q') (20) 

where in the second line we have not written the de- 
pendence on zi, . . . , Zm since by symmetry again it does 
not matter which particles zi, . . . , im we choose to use to 
calculate V^'^- Again, we use the subscript M here to 
indicate that we are considering an M-body interaction. 
Note that this matrix element is necessarily diagonal in 
L by rotational invariancei^ of the interaction V , but is 
not generally diagonal in q. For a given L, if there are 
d different basis vectors \L,q), one must then specify the 
full d by d Hermitian matrix V^'lj. 

At this point, it is useful to point out that we need 
not work with an orthonormal basis \L,q). Even given 
a nonorthonormal (but linearly independent) basis \L,r) 
we can construct an orthonormal g-basis as 

\L,q)=J2[R-'/%r\L,r) (21) 

r 

where the hermitian matrix R is given by 

Rry ^{L,r\L,r') (22) 

Thus, we can freely translate our Hamiltonian into an 
arbitrary \L,r) basis as 

^ \L,r;ii,...,iM)V[llj{L,r';ii,...,iM\ 

ii<...<iM L,r,r' 

(23) 
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where now 



Yr,r 



El 



\r', ,r' 



\R- 



mr[\Vmn) (24) 



The form of Eq. [23] is generally the analogue of the Hal- 
dane Hamiltonian for M-body interactions, and is appli- 
cable for any complete linearly independent basis \L,r) 
(although it may in fact be convenient to choose to work 
with an orthonormal basis). 

While this result seems relatively straightforward, it 
is somewhat nontrivial to find a basis \L,r) (or even to 
compute its dimension). This is precisely the task we will 
undertake in the next section. 



IV. SPACE OF FIXED ANGULAR 
MOMENTUM RELATIVE WAVEFUNCTIONS 

We would like to examine the space of M-particle 
translationally invariant wavefunctions ip with relative 
angular momentum L. In other words, we want to study 
analytic functions in M variables which are homoge- 
neous of degree L and are antisymmetric for the case 
of fermions, or symmetric for the case of bosons. 

Let us start by quoting a result to be proven below 
which applies directly to the bosonic case: The dimen- 
sion of the space of of translationally invariant symmet- 
ric polynomials in M variables which are homogeneous 
of degree L is given by 

dsy,n{L, M) = number of partitions of the integer L 
into pieces no larger than M which 
do not include the integer 1 (25) 

For example, we have (isj,m(7,4) = 2 since there are 
only two partitions of 7 into pieces larger than 1 and 
no larger than 4 (These partitions are obviously 7=3-1-4 
and 7=3-1-2-1-2). A table of dsym{L, M) is given explicitly 
for small values of M and L in table [H In Appendix [C] 
an explicit analytic form is given for this quantity. The 
proof of Eq. [25] will be performed below by construct- 
ing an explicit basis for these symmetric homogeneous 
translationally invariant polynomials. 

To handle the fermion case, it is useful to use the 
result that any antisymmetric polynomial can be writ- 
ten as a Vandermonde determinant times a symmetric 
polynomiaU^ 



'^antisym 

(zi . . . zm) 



i<j 



ZM) 



(26) 

Thus, there is a one to one mapping from the space of 
homogeneous symmetric polynomials of degree L = p in 
M variables to the space of antisymmetric polynomials of 
degree L=p+ {M{M - l))/2 in M variables. It is there- 
fore sufficient to consider only the symmetric case, and 
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ym{L, M), the dimension of the space of 



TABLE I: Table of d. 
homogeneous symmetric translationally invariant polynomi- 
als of degree L in M variables. The formula for these entries 
is Eg. ICll or equivalently Eq. 1251 The dimension of the space 
of homogeneous antisymmetric translationally invariant poly- 
nomials of degree L in AI variables is given in terms of this 
table by danUsym{L, M) = dsym{L - M{M - 1)/2,M) (see 
Eq. I27|) . Note that for an A/-body interaction, at angular 
momentum L, there is a d by d Hermitian pseudopotential 



matrix (d 
fermions) . 



dsym{L,M) for bosons and da-ntisym{L, M) for 



the dimension of antisymmetric homogeneous functions 
of degree L in M variables is given by 

^antisym 

(i,M) =4ym(i-M(M-l)/2,M). (27) 

We now proceed to construct an explicit basis for the 
homogeneous translationally invariant symmetric poly- 
nomials of degree L in M variables. To do this, we begin 
with the well known elementary symmetric polynomials^^ 



..AI 



.M 



(Zl, . . . ,Zm) 



0<ii<i2<... 



(28) 

for m < M and em,M is defined to be zero for m > M. 
The products of these elementary polynomials form a lin- 
early independent basis for all symmetric polynomials in 
M variables. (I.e., the elementary symmetric polynomi- 
als generate the mathematical ring of symmetric polyno- 
mials.) We can write this basis explicitly as 



M 



VXM 



Y[ [ej,Mi^i^ ■ ■ ■ ,^Af)]^ 



(29) 



where A is a vector of nonnegative integers Xj with 
j — 1, . . . ,M. The degree of the basis vector vx^m is 

L = X]jf=i j^j- Thus, we can think of each different A as 
representing a partition of the integer L into pieces with 
each piece no larger than M. For example, the vector 
A = (2, 1, 2, 3) represents the partition (1, 1, 2, 3, 3, 4, 4, 4) 
which corresponds to a homogeneous symmetric polyno- 
mial of degree L = 22. Thus the dimension of the space 
of symmetric homogeneous polynomials of degree L in 
M variables is given by the number of partitions of the 
integer L into pieces no larger than M (denoted Pl,m in 
appendix [C| . 
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The subset of polynomials we are interested in are, in 
addition, translationally invariant. Any such polynomial 
can be written as a function only of the relative coordi- 
nates 

1 ^ 

" " M ^ (30) 

We thus consider as our basis polynomials generated 
by products of em.Aiizi, . . . , zm) Note however that 
since ei^/ = ei^izi, ■ ■ ■ , zm) = Y.f=i '^e have 
ei,Af(2i, • ■ • , zm) — 0. Thus, our new basis is of the form 

M 

wx,M = Yl . . . , zm)]^' (31) 

which for a given A is of degree L = J2jL2 ^^j- Analogous 
to above A now represents a partition of L into pieces no 
larger than M but now none of the pieces are allowed to 
be 1 (I.e, we can take Ai — 0). The size of this basis 
is therefore precisely given by the rule quoted above in 
Eq.[2ni 

The wx^M are constructed to span the space of trans- 
lationally invariant symmetric polynomials. However, we 
must also prove that the all of the wx_m basis elements 
are linearly independent (This is nontrivial since there 
is a linear relation between all of the z's). This proof is 
more involved and is relegated to appendix |DJ 

V. DISCUSSION 

As in the original Haldane^ construction of pseudopo- 
tentials for two-body interactions, one of the main uses 
of our result is to describe any interaction in a maxi- 
mally simplified parameterized form. However, as with 
the Haldane case, examining the form of the pseudopo- 
tential coefficients (and in particular, examining table 
m can have some interesting ramifications in terms of 
what type of simple quantum Hall states might exist. 
To see this connection, let us revisit the Laughlin states 
for a moment. As pointed out by Haldane's original 
work^, the Laughlin state can be thought of as the ex- 
act highest density zero energy state of an interaction 
with only certain nonzero pseudopotentials. For simplic- 
ity here, let us continue to think about quantum Hall 
effect of bosons (i.e., symmetric quantum Hall wavefunc- 
tions). If we consider two-body potentials, examining 
table [J we see that dsym{L, 2) = 1 for even L and is zero 
for odd. Thus there is a single pseudopotential coeffi- 
cient for even values of L, and no coefficient for L odd. 
I.e., to specify a two-body potential we specify the co- 
efficients Vb.2, ^2,2, V4.2, . ■ . (The traditional notation is 
Vo, V2, V4, .... Recall the additional subscript 2 here 
refers to the fact that we are considering a two-body in- 
teraction). Haldane's interesting realization^ is that the 
Laughlin bosonic z/ = 1/2 state is the exact highest den- 
sity zero energy state of a potential that has Vb,2 > 



and all other V's equal to zero. Similarly, if Vo,2 > and 
V2.2 > 0, and all other coefficients are zero, the highest 
density zero energy state is the Laughlin bosonic v = 1/A 
state. These results are, of course, well known. We em- 
phasize that to obtain these Laughlin states we are sim- 
ply forbidding any two particles from having a relative 
angular momentum below some cutoff. 

Let us now consider three-body interactions for bosons. 
Examining our table, we see that dsymiL,3) = 1 for 
L = 0,2,3,4,5,7. Thus, for each of these values of L, 
there is a single pseudopotential coefficient (for L — 1 
there is no state, since dsym{^,3) = and therefore no 
coefficient). In analogy with Haldane's original paper, 
one might naturally ask questions like what is the high- 
est density zero energy state of a system with V0.3 > 
and all other pseudopotential coefficients equal to zero. 
In other words, what happens if we forbid any three 
particles from having zero relative angular momentum. 
(Note that this Vb,3 pseudopotential coefficient also cor- 
responds to a three-body point-contact interaction, as 
demonstrated in appendix IE)) . It is well known that the 
highest density zero energy state of such an interaction 
is the bosonic Moore-Read Pfaffiar4ii^ at filling fraction 
1^ = 1. 

We could similarly ask what happens if we make Vb,3 > 
and ^2,3 > 0, and all other coefficients zero. Thus 
we are asking what happens if we forbid three particles 
from having relative angular momentum or 2 (there 
is no state where three particles have relative angular 
momentum 1). The highest density zero energy state 
of this interaction, known as the "Gaffnian", occurs at 
filling fraction = 2/3 for bosons and has been studied in 
depth in Ref.|7(. (In appendix lEl a form of this interaction 
is given as derivatives of delta functions). Similarly, if 
M3.3, ^2,3, and V3^3 are positive and all other coefficients 
are zero the highest energy zero energy wavefunction of 
such an interaction is known as the Haffnian and has 
been studied in depth in Ref. 0. In each case here, we 
are simply forbidding three particles from having relative 
angular momentum below a certain cutoff and we can 
consider successively increasing this cutoff. 

One can also generalize this construction to other val- 
ues of M. For example, the Read-Rezayi stated at filling 
fraction 1/ = {M — l)/2 for bosons, is the highest density 
zero energy state of a Hamiltonian that consists only of 
a positive Vo,m coefficient and all other coefficients are 
zero. The generalization of this construction — forbid- 
ding clusters of M particles to have angular momenum 
below some cutoff L — is studied by the current authors 
in Ref. B 

When dsym{L, M) > 1 there 

a- dsym dimensional pseu- 
dopotential matrix V^'\j. For simplicity, let us assume 
that we are working in a basis where the wavefunctions 
\L,q) are orthogonal states (recall that each \L,q) is an 
M-particle relative wavefunction). If we were to choose 

the V^'lf matrix to be positive definite (all eigenvalues 
strictly positive), it would forbid (give positive energy 
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to) any cluster of M particles with relative angular mo- 
mentum L. This is the situation that is considered in 
Ref. i. 

However, we could also consider more interesting sit- 
uations. If we choose V^'Xf to have some eigenvectors 
with positive eigenvalue and some eigenvectors with zero 
eigenvalue, we allow only certain clusters of M parti- 
cles with angular momentum L while disallowing (giv- 
ing positive energy to) others. To be more specific, any 
wavefunction describing a cluster of Af -particles that lies 
completely in the space of zero eigenvalue eigenvectors 
(the "null space") is a zero energy cluster, whereas any 
cluster that has components outside of this null-space has 
positive energy. (An explicit example of this is given in 
Appendix [E]). 

It is particularly interesting to consider the situation 
where the matrix V^'^ has a single eigenvector |(/)) with 
zero eigenvalue (and all other eigenvalues are positive). 
What this means is that in the dsymiL, M) dimensional 
space of translationally invariant homogeneous polyno- 
mials of degree L in M- variables, there is only a single 
(translationally invariant) polynomial (j){zi, . . . , zm) that 
is a zero energy M-particle wavefunction, and all other 
polynomials are positive energy wavefunctions. We may 
then ask whether a quantum Hall wavefunction exists 
with the property that all clusters of M particles are in 
the relative wavefunction (f>. If such a wavefunction does 
exist, it would certainly be a zero energy state of our 
chosen Hamiltonian. In some cases it may also be the 
highest density zero energy state, i.e., the ground state 
of the Hamiltonian. In future work, we will show a num- 
ber of examples where this is indeed the case. 

Such a general construction could be quite significant, 
as it allows us to construct a different quantum Hall wave- 
function for each (translationally invariant) cluster func- 
tion (p that we might define. In cases where the dimen- 
sion dsym{L, M) — 1, this is no freedom at all - there is 
only a single possible cluster function 0. However, when 
dsym > 1, we have a continuum of possible choices — and 
therefore can define a continuum of possible Hamiltoni- 
ans with a continuum of ground states wavefunctions. 

Most interestingly, by judiciously choosing the clus- 
ter function we may be able to find Hamiltonians 
whose ground state are quantum Hall wavefunctions cor- 
responding to known unitary conformal field theories 
(For discussion of the connection between quantum hall 
and conformal field theory see, Refs. HIslflTh. Without 
proof we believe such constructions are generically pos- 
sible, and we will show several examples of this in future 
work. The general belief is that any wavefunction corre- 
sponding to a unitary conformal field theory corresponds 
to a stable phase of matter. 

To summarize our work, the main purpose of this pa- 
per was to formulate Haldane pseudopotentials for multi- 
particle interactions. The main technical result of this 
paper is the enumeration of the possible relative wave- 
functions for clusters of M-particles. Since we only con- 
sider interactions that are rotationally invariant, we cat- 



egorize the possible M-particle translationally invariant 
cluster wavefunctions by their total angular momentum. 
For bosons, we find dsym{L, M) basis states whereas for 
fermions we calculate dantisym{L , M) . We discuss how 
a rotationally invariant A/-body interaction can be de- 
composed into pseudopotential coefficient matrices V^'X/- 
For a given L and M the matrix is d by d (with d being 
dsyrn for bosons and dantisym for fermions. We further 
describe a relatively simple linearly independent (but not 
orthonormal) basis in which such an interaction matrix 
could be expressed. 

For increased clarity, we would like to show the some 
simple examples of using this pseudopotential formalism. 
In particular, we will give details for bosonic systems with 
three-body potentials. However, since these examples are 
somewhat lengthy we have put them in Appendix [El We 
encourage the reader to work through the calculations 
given there to cement the general ideas of this paper. A 
brief summary of the appendix is as follows: 

• We write down the wa.3 basis for three-cluster 
wavefunctions (as discussed in section IIVI above) 

• Since this basis is not orthonormal, we show how 
to orthonormalize (as discussed in section WU^ 

• We consider several simple toy model three-body 
interactions and write down their pseudopotential 
coefficients in this basis (as discussed in section lHl)) 

• We show how to calculate pseudopotential coeffi- 
cients for higher Landau levels (analogous to that 
discussed in section HH 

The explicit demonstration of these steps in the appendix 
should make concrete the abstract ideas expressed above. 
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APPENDIX A: SPHERE AND TORUS 
GEOMETRIES 

On the spherical geometry^, angular momentum is a 
good quantum number, and any particle in the LLL has 
angular momentum S — with the total number 

of flux quanta through the sphere. The maximum total 
angular momentum of M particles is MS (if they are 
bosons). What we called "relative angular momentum" 
L on the plane corresponds to total angular momentum 
MS — L on the sphere. For fermions, since the minimum 
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relative angular momentum is M{M — l)/2, the max- 
imum total angular momentum is MS — M{M — l)/2. 
The direction of the "total" angular momentum plays the 
role of center of mass coordinate on the plane. Interac- 
tions that are rotationally invariant, again, will not mix 
states of differing "relative" angular momentum. 

A way to think about the translation from plane to 
sphere is to realize that any M-particle cluster wavefunc- 
tion is completely determined by its behavior as all of the 
particles come close together — the "analyticity" of the 
wavefunction implies that we know the complete wave- 
function over the entire system if we know its behavior 
for short distances. Thus, if we look at M particles in 
some relative wavefunction very close together, they do 
not "see" the curvature of the sphere — locally we can 
think of the particles as living on a plane, and so our un- 
derstanding of the possible wavefunctions on the plane 
can tell us all the possible wavefunctions on the sphere. 
Once we specify the analytic form of the wavefunction for 
the particles very close together, this uniquely specifies 
the wavefunction over the entire sphere, thus generating 
an appropriate basis for "relative" wavefunctions on the 
sphere. In other words, for each of our relative wave- 
function states on the plane, there is exactly one relative 
wavefunction on the sphere with the same analytic form 
at short distances. 

In fact, it turns out that there is an extremely simple 
translation from the plane to sphere. If we point our 
center of mass coordinate towards the north pole, there 
is then a simple mapping of the planar wavefunctions 
(with center of mass at position zero) to the sphere, by 
stereographic projection^. The only slight complication 
is that the usual gaussian factors on the plane become 
modified 

p-|^lV4^ ]. (-An 

(1+ |z|2/4il'2)l+S ^ > 

where R is the radius of the sphere. 

On the torus, the story is somewhat more complicated, 
but in principle similar. First of all, on the torus angu- 
lar momentum is not conserved (for arbitrary unit cell). 
Nonetheless, we can define a linearly independent basis 
for relative wavefunctions. As with the above discussion 
for the sphere when we bring the particles close together, 
the systems looks locally just like an infinite plane. For 
each state in our linearly independent basis for the plane 
there is exactly one state on the torus that has the same 
short distance limit. Thus, in principle, we can convert 
our basis on the plane into a basis on the torus. 



APPENDIX B: REDUCIBLE INTERACTIONS 

Given an M-body interaction ^(ri, . . . ,tm) it is fre- 
quently convenient to separate out any part of an interac- 
tion that is "reducible" to interactions between any lower 
number M' < M of particles. For example, a three-body 



interaction of the form 

y(ri,r2,r3) = V2(ri-r2)+V2(r2-r3) + V^2(ri-r3) (Bl) 

is fully reducible to two-body components. 

For an arbitrary M-body interaction, it is easy to find 
the M'-body reducible component (with M' < M) by 
taking appropriate limits where M — M' of the particles 
are moved away to infinity. For example, for a three-body 
interaction, the two-body part is given by 

^2(ri-r2)= lim l/(ri,r2,r3) (B2) 

We then have the irreducible three-body interaction given 
by 

Virriri,V2,rz) = V{yi,T2,r:i) (B3) 
- [V2{VI - V2) + V2{V2 - r3) + V2{yi - r3)] 

We note that it is possible (although unnecessarily com- 
plicated) to describe any two-body interaction as an 
unreduced three-body interaction. However, typically it 
is simpler (and clearer) to fully reduce three-body inter- 
actions and treat their two-body piece separately. For 
most realistic potentials, the two-body piece will be very 
large compared to the reduced three-body piece (which 
in turn will be large compared to the reduced four-body 
piece, and so forth). However, by fine-tuning, in cer- 
tain situations it can be arranged that the pure two-body 
piece vanishes-. 



APPENDIX C: ANALTYIC EXPRESSION FOR 
PARTITIONS 

We claim that Eq. [25] can be reexpressed as 

dsy,n{L, M) = Pl,m - Pl-i.m (CI) 

where Pl,m is the number of partitions of the integer L 
into pieces no larger than M. (Note also that Pl.m is 
equivalently the number of ways to partition the integer 
L into no more than M pieces). 

To see the equivalence between Eq. ICll and Eq. [21] we 
simply note that given any partition of L into pieces no 
larger than M which includes the integer 1, we can drop 
this 1 and (uniquely) generate a partition of L — 1 into 
pieces no larger than M . A generating function for Pl,m 
is given analytically as 

M 00 

= n - E 'i'p^M (C2) 

m=l ^ L=l 

so that 
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APPENDIX D: PROOF OF LINEAR 
INDEPENDENCE 

Linear independence of the wx^m basis can be refor- 
mulated as follows: Given a polynomial P in the space 
spanned by the wx^m (i-e, given a translationally invari- 
ant symmetric polynomial in M variables) this polyno- 
mial can be written as 

P = ^ OA Wx,M (Dl) 
A 

where Ai = in all terms of the sum (since the wx.ai 
basis spans the space) . Proof of linear independence of all 
of the wx,M basis vectors is equivalent to the statement 
that each such polynomial corresponds to a unique set of 
coefficients ax- 
To prove this linear independence, we make use of 
translational invariance again and shift the center of mass 
coordinate of all M of the z^'s to zero. Thus, we can 
rewrite Eq. [31] as 

M 

Wx,M = Y\_ [ej,M(2:i, • ■ ■ , ZM-1, z)]^' (D2) 

where 

M-l 

z=-J2^j (D3) 

Here, wx.m must be a symmetric polynomial in the re- 
maining M — 1 variables Zj. Indeed, it is easy to check 
that 

e7n,M{zi, ■ ■ ■ ,zm-i,z) = emM-i{zi, . . . ,zm-i) 
(zi, . . . ,zm-i) ei,M„i(zi, . . . ,zm-i) 
— e^.M-i — e.m-i,M-i ei,M-i (D4) 

where again e^.M-i is defined to be zero for m > M . In 
particular note that eM,M-i = and 

eM,Af(2:i, • ■ • , ZM~i, z) = eM-i,M-iei,Af-i (D5) 

By rewriting our basis vectors wx,m in terms of the 
em>/-i we will be able to use the fact that the basis 
vectors generated by these em,M-^ (i-*^-' ^^a.m-i ba- 
sis) are linearly independent. 

We proceed as follows: Given the polynomial P written 
in terms of the wx basis vectors as in Eq. IDl) we rewrite 
each wx as sums of products of em,^/-! using Eq. ID4I 
I.e., we are rewriting everything in terms of the wa.az-i 
basis. We then group terms such that P — Pq + Pq 
where Pq includes terms with no factors of e\,M~\ f-nd Pq 
includes terms with at least one factor of e\^M-\- Writing 
Pq in the f A,Af-i basis we have 

A(0) 



where since Pq has no factors of e\,M-\^ the vectors A*^^ 
must have A^^"* = 0, and we must also have Am'' = 
for m > M since Cm.Af-i = for m > M. We then 
construct, 

A(0) 

with the basis wx,m defined as in Eq. ID2I Examining 
Ea. lD41 we see that the e„_M-i terms generate precisely 
the polynomial Pq within the expansion of Qq. The other 
terms in the expansion of Qq all contain ci^m-i- Note 
that unlike Pq, we have constructed Qq to be in the space 
spanned by wx,m (i.e., it is a translationally invariant 
polynomial in M variables). Thus, since P is in the space 
spanned by wx,m then P—Qo is also in the space spanned 
by Wx.m- The difference P — Qo is constructed so as to 
have ei M-i in a-U of its terms. Examining Eq. ID4I and 
Eq. ID5I we see that if a polynomial (P — Qq in this case) 
is in the space spanned by wx.m and it has ei,M-i in all 
of its terms, it must also contain a factor of cm-laz-i 
(i.e., this factor could only have come from a factor of 
CAf.Af in Wx.m) or it must be zero. Thus, we can write 

P = Qq - ei^MCAf-l.Af-lPl (D8) 
= (3o + eAf,Af(zi,...,ZAf-i,z)Pi (D9) 

and now Pi is a polynomial of degree less than that of 
P and is in the space spanned by wx,m- We then iterate 
this procedure, to rewrite Pi as a sum of two terms, Qi 
containing no factors of ei jv/-i and another term eM,A/P2 
with P2 of lower degree still, and so forth. Eventually for 
some Oimaxi we must obtain Pama^ ^ Qama^ = 0; which 
then terminates the procedure. Thus, we can successively 
decompose into a finite sum 

[eA/,Af(^l, • • • jZM-I^z)]" Qa (DIO) 

a=0 

Where at each level of the iteration, we have 

A<°) 

where the vectors A*-"-* must have A^"'' = A^"'' = 0. We 
realize that Eqs. IDlOl and lDlll can be recast in the form 
of Eq. IDll where the coefficients ax are determined by 
ax — o,^^(L) when Xm — a and a[,"' — \k for all k ^ M 

(and we always have A^"'* = Ai = 0, and Am = A,„ = 
for m > M). 

Our construction here is unique in the sense that given 
any polynomial P in the space spanned by wx,m we find 
a unique set of coefficients ax in Eq. IDll This then com- 
pletes our proof of the linear independence of the basis 

WX.M- 
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APPENDIX E: EXAMPLES: BOSONS WITH 
THREE-BODY INTERACTIONS 

To elucidate our results, we will explicitly consider 
bosons interating via three-body interactions. As dis- 
cussed at length above, we will be concerned with three- 
body wavefunctions. In calculating any matrix element 
of these wavefunctions, we mean 

(El) 

where O is any operator and the integration measure is 

dfi = dzi dzi dz2 dz2 dz^ dz^ (E2) 

Recall that our full wavefunctions include both the trans- 
lationally invariant relative wavefunction ip as well as and 
the center of mass coordinate and the Gaussian 

factors (See Eq. [HI) 

(E3) 

Assuming that O is translationally invariant, it only cou- 
ples to Tp^'^'K Thus, we can choose an orthogonal basis 
for ijj^'^'^ and we will have the matrix element be zero if 
is orthogonal to ip^^ . If we choose *^ = V'?*^ 
then, since the operator is translationally invariant, the 
matrix element should be independent of our different 
possible choice of -0'^*^, so long as the different possi- 
ble choices are normalized the same way. For performing 
calculations, it is then acceptable to choose, for exam- 
ple, ip^'^^ = tp2^ = 1, and we can then check that our 
results are independent of this choice. This also implies 
a normalization choice for our matrix elements. We will 
work with this choice throughout this appendix. 

For the relative wavefunctions, we will work in the 
wx,M basis discussed above. We write out the first few 
three-cluster wavefunctions here explicitlyi^ 

\L = 0) = i«(o,o,o),3 = 1 (E4) 

\L = 2) = W(o, 1,0), 3 = 62,3(^1,22,53) 

\L = 3) = W(o, 0,1), 3 = 63,3(21,22,23) 

|i = 4) = W(o,2,0),3 = [62,3(^1, 22, 23)]^ 

\L = 5) = W(o, 1,1), 3 = 62,3(^1,^2,23)63,3(51,22,23) 

|L = 6,a) = ^(o, 3,0), 3 = [62,2(21,22,23)]^ 

\L^6,b) ^ u'(o,o,2),3 = [62,3(21, 22, 23)]^ 

\L = 7) = ^(o, 2,1), 3 = [62,2(21, 22, 23)1^63,3(21, 22, 23) 



and recall the definition, Eq[2Dl of Zi in terms of 2^. Note 
that none of these basis states are normalized. Further, 
note that since dsymiL = 6,3) = 2 there are two basis 
states for L = 6. The two states, as discussed above, are 
linearly independent, but not orthogonal (or normalized). 



We remind the reader that a wavefunction with relative 
angular momentum L vanishes as L powers when all three 
particles come to the same point. 

We note in passing that for the particular case of three 
body wavefunctions, another basis has been constructed 
that is already orthonormalizedi^. However, this con- 
struction does not generalize easily to more than three 
particles, so we choose to work with the above basis to 
show the more general method. 

We start by calculating the normalization of these basis 
states, by setting the operator O — 1. The integrations 
over 2i's can then be performed straightforwardly (Math- 
ematica makes this trivial) to give the normalizations 

{L = 0\L = 0) = 2^n^ 

{L = 2\L = 2) = 2-V^ 

{L = 3\L = 3) = 2^7rV3^ 

(L = 4|L = 4) = 2^n^ 

{L = 5\L^5) = 2i2^V3' 

{L^7\L = 7) = 5*2^^7T^/3^ 



(E5) 



and for L — 6 the full overlap matrix (See Eq. 

by 



IS given 



{b\a) m 



32*2" 



-212/3 
212/3 11*213/3"* 



(E6) 

Example 1: As an example we now consider the toy 
model three-body point contact interaction 

F(ri,r2,r3) = V d{ri - r2)(5(r2 - rg) (E7) 

Note that this three-body interaction is fully reduced in 
the sense of Appendix |B] Allowing the delta functions to 
act, the interaction matrix element then reduces to 

(^'i|-^^|*2) = vJdzdz*[Mz,z,z)]*Mz,z,z)e~^^'^"/^ 

(E8) 

From our basis states, only the state |L = 0) = 1 does not 
vanish when all three particles come to the same position, 
so the only nonzero matrix element is trivially calculated 
to be 



{L = 0\V\L = 0) = (27r/3)f 



(E9) 



Thus, the only nonzero pseudopotential coefficient is Vo,3. 
In our non-normalized basis (See Eq. I24p we obtain 



^0,3 = 



V 



{L^Q\V\L = 0) _ 
|(L = 0|L = 0)|2 ~ 3(27r) 



(ElO) 



Alternately, we could construct a normalized basis (See 
Eq.ini) 



\L = 0)no 



\L = 0)/V8^ 



(Ell) 



in terms of which we would have (See [20]) a pseudopo- 
tential coefficient 



Vo,; 



No 



,{L = 0\V\L^O) 



V 



Nc 



3(2^)2 



(E12) 
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As mentioned above, the highest density zero energy 
state of this interaction is precisely the Moore-Read Pfaf- 
fian state^ii^. 

Example 2: Let us now consider a more complicated 
interaction. Consider 

V{yuY2, ra) = V V\[ 5{v^ - ra) 5{y2 ~ rg) ] (E13) 

where the subscript 1 on means we are differentiat- 
ing with respect to the position ri. Derivative of delta 
function interactions, such as this one, are quite useful 
since arbitrary interactions can be built up as a series of 
successive derivatives of delta functions (We can think 
of such a series as a Taylor series expansion in fourier 
space). 

To handle this interaction, we first integrate by parts 
then let one delta function act to give 



*i|V"|*2) = V / dridr2(5(ri - ra) 



V\ [ll^l{zi, Z2, Z2)]*1p2{zi, Z2, Z2)e 

Using the fact that 

the second line of Eq. IE14I then becomes 

(42)[(9,. -zi/2)2[V;i(zi,Z2,^2r] X 
[id,,-zl/2fMzi,Z2,Z2)] e~l"^l'/' 



/2 



(E14) 



(E15) 



(E16) 



Once the last delta function acts, all three of the particles 
are put at the same position. As we mentioned above, 
the wavefunction for a state with angular momentum L 
vanishes as L powers when all the particles come to the 
same position. Now, since we have up to two derivatives, 
we see that both L = and L = 2 can have nonzero 
matrix elements, but no higher L wavefunction can. The 
L = matrix element is messy to calculate (although it 
is easy on Mathematica) and gives 



{L = 0\V\L = 0) = {4^tt/27)V 



(E17) 



The L = 2 matrix element, on the other hand, is actually 
easy to calculate. Here, we only need to keep the terms 
where all of the derivatives act on the wavefunctions (the 
z/2 and z* /2 do not contribute). With our wavefunction 
being 62,3, we note that 



9^162,3(^1,^2,23) = 
9^162,3(^1,^2,23) = 

so that 



di^e2,3{zi,Z2,Z3) = iz2 + Z3) 

-2/3 (E18) 



(-2/3)|i = 0) 



(E19) 



Thus the matrix element is given by (—2/3)^ times the 
value of Eq. IE9I times the prefactor of 4^ from Eq. IE15I 



Again, we could write this in either the normalized, or 
unnormalized basis to give the pseudopotential coefficient 
^2,3. 

Example 3: Let us now try a much more complicated 
interaction. Consider 

l/(ri, r2, r3) = F [^(ri - r2) (5(r2 - r3) ] (E21) 

The reason we choose this is because it has nontrivial 
matrix elements of L = 6, which is the first "interesting" 
case where there are two states at the same L. As above, 
we begin by integrating by parts. As in Example 2 above, 
calculation of matrix elements for L < 6 is quite messy, 
but for L = 6 things simplify quite a bit. Analgous to 
the L — 2 case for Example 2 above, the only term that 
does not vanish when the delta function acts is the one 
where all of the derivatives have been applied. Noting 
that from Eq. IE4I we have \L = 6, a) = [62,3]^ and \L = 
6, b) — [63.3]^ we then have 



/3 



^ <[e2,3] 



6! 



[9.^62,3]^ = -80/3 



2!2!2! 

^[9.^63,3]^ 



320/81 



Where we have used Eq. IE18I above, as well as the anal- 
ogous easily calculated 9^^63,3 = 4/9. With these facts, 
to determine the matrix elements of V we realize that 
we have just the same integral as in Eq. IE9I above times 
these factors (and a prefactor of 4^ coming from Eq. IE15[) 



{a\V\a) {a\V\b) 
{b\V\a) {b\V\b) 



{4^){27t/3)V 



o? a/3 
a/3 0^ 



(E22) 



Had we chosen to consider instead 

T/(ri,r2,r3) = VV\V\ [ 5(ri - r2) (5(r2 - r3)] (E23) 



(L = 2\V\L = 2) = (42)(87r/27)l/ 



(E20) 



we would have obtained the same form of Eq. IE22I but 
with 

a = 93^93^62,3]' = 28/3 (E24) 
/3 = 93^ 93je3,3]' = 104/81 (E25) 

Or analogously we might have chosen the interaction 

y(ri,r2,r3) = Y V\v\v\\b{v^ - X2)&{v2 - v^)] (E26) 

to obtain the form of Eq. IE22I but with 

a ^ dl^dl^dl^{e2,3f = -%lZ (E27) 
/3 ^ 9f^9f,92je3,3]' = 176/81 (E28) 

and so forth (obviously there are many more possibilities 
in the same spirit). 

It is interesting to note that independent of which of 
these interactions we choose, the form of Eq. [E22l is rank 
one — meaning that it leaves some three-cluster wave- 
functions with i = 6 at zero energy as discussed in sec- 
tion |Vl If we want to give positive energy to all three- 
cluster relative states with L = 6, we can do so by adding 



11 



together two of these interactions, for example adding the 
interaction in Ea. lE21l to the interaction in Eq. lE23l (and 
we should choose V positive in both cases). 

Given the result in Eq. IE22I and given the normaliza- 
tion matrix Eq. lE6[ we can plug into Eq.[24]to obtain the 
pseudopotential matrix Vq'^ in this non-orthonormal ba- 
sis. Or, we could construct an orthonormal basis using 
Eq. 1211 and then construct the pseudopotential matrix 
V(3,3 in this basis as in Eq. [20l 

Example 4: Let us now consider bosons in the first 
excited Landau level interacting with a three-body po- 
tential. Note that this is a bit of an artificial problem 
since most interesting boson problems would be lowest 
Landau level. Nonetheless, it is a well defined question 
of what the pseudopotentials would be for bosons in any 
given Landau level. Analogous to Eq. [Til we can raise 
any three-particle wavefunction into the first excited Lan- 
dau level by applying ajojag to the wavefunction, with 
the raising operator given by Eq. [T^ To raise to the ri*'* 
excited Landau level, we would use (a^)"'- 

We can thus raise the entire basis given in Eq. IE4I to 
create a basis in the first excited Landau level, which we 
write as 

\L = 0)iLL = a\alal\L = 0) ^ z^z^z; /2V2 
\L = 2)iLL = a\alal\L ^ 2) = . . . 

: (E29) 

with the unraised states on the right defined as in Eq. lE4l 
We have written out the case of L = explicitly (with- 



out the Gaussian factor), but the case of L = 2 is a 
rather long expression. It is worth noting, however, that 
very generally, the raised function |i)iLL does not van- 
ish when all three particles come to the same position for 
i < 3 (the number 3 occurs here because there are 3 par- 
ticles, therefore 3 raising operators, therefore a maximum 
of 3 derivatives). 

It is convenient that the normalizations of these states 
are unchanged 

ill{L^O\L^O)ill = {L = 0\L = 0) (E30) 
iLi(i-2|i = 2)iLi = (L = 2|L = 2) 



since [a\a] — 1, and LLL states plays the role of the 
vacuum. 

Let us now consider the simple three point delta func- 
tion interaction of Eq. IE7I Allowing the delta functions 
to act we obtain the matrix elements as a single remain- 
ing integral as in Eq. IE8I It is trivial to analytically 
obtain the result for L — and L = 3 (we leave this 
as an exercise for the reader). For L — 2, however, the 
calculation is done with the aid of Mathcmatica. We end 
up with 

^ll{L^O\V\L^O)ill - (47r/27)f (E31) 
ill{L^2\V\L^2)ill - {8n/9)V 

ill{L = S\V\L = 3)ill = (2567r/243)\/ 
and matrix elements for higher L all vanish. 
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